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^ , ABSTRACT 
Pk 

Aims. We investigate the spectral characteristics of 51 candidate High Frequency Peakers (HFPs), from the "bright" HFP 
sample, in order to determine the nature of each object, and to obtain a smaller sample of genuine young radio sources. 
Methods. Simultaneous multi-frequency VLA observations carried out at various epochs have been used to detect flux density 
and spectral shape variability in order to pinpoint contaminant objects, since young radio sources are not expected to be 
significantly variable on such a short time-scale. 

Results. From the analysis of the spectral variability we find 13 contaminant objects, 11 quasars, 1 BL Lac, and 1 unidentified 
object, which we have rejected from the sample of candidate young radio sources. The ~6 years elapsed between the first and 
latest observing run are not enough to detect any substantial evolution of the overall spectrum of genuine, non variable, young 
radio sources. If we also consider the pc-scale information, we find that the total radio spectrum we observe is the result of the 
, superposition of the spectra of different regions (lobes, hot-spots, core, jets), instead of a single homogeneous radio component. 
■ This indicates that the radio source structure plays a relevant role in determining the spectral shape also in the rather common 
^ / ' case in which the morphology appears unresolved even on high-resolution scales. 
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1. Introduction H9951 Readhead et al. [T996l Snellen et al.[2000). 

5-j However, it has also been claimed (Alexander 2000; 

The origin and evolution of radio emission is one of Marecki et al. HDHl that a fraction of young and compact 

the greatest challenges in the study of Active Galactic radio sourceg may die in an early gtage before becoming 

Nuclei. The onset of radio activity is often thought to be j e gcale ob j ects 

linked to merger or accretion events in the host galaxy, The early gtages ^ ^ evolution scheme (typical ages < 

which provide enough fuel to feed the central AGN. As a 1q4 yeargj Polatidis & Conway M M urgia H003D are 

consequence, at least in its early stages, the radio emission represented by the popu i ation of Compact Symmetric 

evolves in a quite dense and possibly inhomogeneous Qbjects (CSQs) They &re & scaled _ down version of thc 

ambient medium which can influence its growth. larg6j powerM radio sources ( cor e-jets-lobes; Wilkinson 

The ideal targets to understand such phenomena are the ^ aJ wM & ^ defined pe&k ^ thdr radio 

young radio sources, whose radio lobes still reside within spectrum at freque ncies between ~ 100 MHz to a few 

the innermost region of the host galaxy. GHz 

Among them, the smallest (i.e. youngest) objects are j n framework, the anti-correlation found between 

the most suitable to investigate the role played by the thfi pe&k frequency and the SQurce size (cmea 

host galaxy Interstellar Medium (ISM) on the evolu- & Baum ^ implies that the youngest sources must 

tion/growth of the radio source. be gought among thoge ^ gpectral peakg occurring 

The evolutionary stage of the powerful radio sources &t frequencies higher than 5 GHZi and termed " H igh 

is related to their linear sizes. Following self-similar Frequency Peakers » (HF Ps, Dallacasa WM- 

evolution models, the most compact sources will evolve The selection of a statistica n y complete sample made of 

into the extended radio source population (Fanti et al. genuine yQung radiQ sourceg Qnly ig essential to study in 
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detail the physical nature of this class of objects. 
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Date 


Conf 


Obs Time 


code 


Sep 12 2003 


AnB 


240 


a 


C 1 -I i \ -"www) 

Sep 13 2003 


AnB 


240 


b 


Sep 14 2003 


AnB 


240 


c 


Sep 15 2003 


AnB 


620 


d 


Jan 22 2004 


BnC 


120 


e 


Jan 26 2004 


BnC 


180 


f 


Jan 28 2004 


BnC 


240 


g 


Jan 30 2004 


BnC 


120 


h 


Mar 21 2004 


CnD 


120 


i 



Table 1. VLA observations and configurations. The total 
observing time (Column 3) is inclusive of the scans on the 
HFPs from the "faint" sample. 



A sample of candidate HFPs (the "bright" sample; 
Dallacasa et al. 2000) was constructed on the basis of 
both the shape of the simultaneous radio spectrum and 
the turnover frequency. However, such selection tools can 
introduce a contamination by beamed radio sources. For 
example, blazar objects, although usually characterized 
by a flat spectrum, can occasionally show a convex 
spectrum when their radio emission is dominated by a 
flaring, self-absorbed component (e.g. Torniainen et al. 
12005! . 

During the majority of their lifetime blazars and young 
radio sources display very different characteristics: the 
former posses significant flux-density and spectral vari- 
ability, and the emission, often polarized, has a Core- Jet 
morphology on various scales (from tens of kpc down to 
pc-scale). The latter have no variability and the radio 
emission, almost completely unpolarized at least at low 
frequencies, has a "Double/Triple" structure (Orienti et 
al.HH]). 

In this paper we present the results of simultaneous multi- 
frequency VLA observations, made at several epochs, 
in addition to those already published by Dallacasa et 
al. ([2000]) and Tinti et al. (f2U05|) . These new observa- 
tions were carried out to further identify blazars with 
variability on time-scales longer than the time elapsed 
between the previous two epochs (Tinti et al. 120051) . and 
to study the evolution of the radio spectrum of those 
HFPs already confirmed as genuine young radio sources 
(Orienti et al.HOUS]). 

In general, they better characterize the properties of the 
sources of the "bright" HFP sample. 



2. Multi-frequency VLA observations and data 
reduction 

Multi- frequency VLA observations of 51 of the 55 can- 
didate HFPs were carried out in different runs from 
September 2003 to March 2004, using filler time. About 



25 objects were observed twice, in order to identify 
short-period variability. 

The observing bandwidth was chosen to be 50 MHz per 
IF. Separate analysis for each IF in L, C and X bands 
was carried out to improve the spectral coverage of the 
data, as done in previous works (Dallacasa et al. [2000; 
Tinti et al. 12005]) . 

We obtained the flux density measurements in L band 
(IFs at 1.465 and 1.665 GHz), C band (4.565 and 4.935 
GHz), X band (8.085 and 8.465 GHz), U band (14.940 
GHz), K band (22.460 GHz), and for a few datasets in Q 
band (43.340 GHz). 

Each source was typically observed for 50 seconds at each 
band, cycling through frequencies. Therefore, the flux 
density measurements can be considered simultaneous. 
For each observing run and at each frequency, about 3 
min were spent on the primary flux density calibrator, 
either 3C286 or 3C48. Secondary calibrators, chosen to 
minimize the telescope slewing time, were observed for 
1.5 min at each frequency every ~20 min. This implies 
that the observations do not have astrometric accuracy. 
Accurate positions of the sources can be found in the 
JVAS catalogue (Patnaik et al. 1992] Browne et al. 
IT998] Wilkinson et al. IT998]) . Information on the date 
and duration of the observing runs is summarized in 
Table [T] The total observing time also takes into account 
scans on candidate HFPs from the "faint" HFP sample 
(Stanghellini et al., in preparation), in order to optimize 
the observing schedule. 

The data reduction was carried out following the standard 
procedures for the VLA, implemented in the NRAO AIPS 
package. Images for each frequency were produced. In 
order to obtain accurate flux density measurements in 
the L band, it was necessary to image several confusing 
sources falling within the primary beam, and often 
accounting for most of the flux density within the field of 
view. 

During some observing runs, strong RFI affected the 1.665 
GHz data, precluding the flux density measurements for 
a few sources at such frequency. 

The final images were produced after a few phase-only 
self-calibration iterations, and source parameters were 
measured by means of the task JMFIT, which performs 
a Gaussian fit. The integrated flux density was checked 
with TVSTAT. The flux density measurements at each 
frequency and epoch are reported in Tabled] 
Apart from a few sources that were already known to 
possess an extended emission (Tinti et al. I2005[) . the 
majority of HFPs are unresolved with the VLA, even in 
the Q band. 

The r.m.s. noise level on the image plane is not relevant 
for bright radio sources as our targets. In this case, the 
main uncertainty comes from the amplitude calibration 
errors, which is within (la) 3% in L, C and X band, 5% 
in U band, and 10% in K and Q band. 
With these new- multi-frequency observing runs, 48 out of 
the 55 sources from the bright HFP sample have at least 
three epochs of VLA data. The comparison between each 
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Fig. 1. Radio spectra of the 51 candidate HFPs observed with the VLA during the observing runs presented in this 
paper [S v (Jy) vs v (GHz)]: crosses, plus signs and diamonds refer to the first (1998-1999; Dallacasa et al. 2000), 
second (2002; Tinti et al. 12005 j) . and our third epoch (2003-2004; dates are reported in Table[T]) of simultaneous VLA 
observations, while the dashed alternated by three dots, the dotted and the solid lines indicate the corresponding fits. 
When a fourth epoch was available (2004; dates are reported in Table [I]), it is represented by triangles and a dash-dot 
line. 
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epoch allows us to determine the variability properties of 
each single source. 



3. Results 

The lack of any spectral variability is a key element for 
the selection of genuine young radio sources. The new 
epochs of simultaneous multi-frequency VLA observa- 
tions, carried out 5-6 years after the first observing run 
(Dallacasa et al. 2000), provide further information on 
the spectral shape and variability of the candidate HFPs, 
allowing us to better discriminate between genuine young 
radio sources and contaminant objects. 



3.1. Spectral shape 

The shape of the radio spectrum is one of the key 
identifying characteristics of young radio sources. The 
overall shape of their radio spectra is convex with a peak 
at high frequencies that is likely due to synchrotron self- 
absorption within the small radio emitting region (Snellen 
et al. I2000j) . although free- free absorption (Bicknell et al. 
119971 Kameno et al. I2000p may play a role. 
On the other hand, blazar objects, which usually show 
flat spectra, sometimes can display a convex spectrum as 
a consequence of an increment of the flux density, when 
a single, homogeneous and boosted component at the 
jet base temporarily dominates the radio emission. As 
the component adiabatically expands, the flux density 
decreases, leading the spectrum back to a flat shape. 
Following the approach from Dallacasa et al. (j2000j) and 
Tinti et al. (200.5), we fit the simultaneous radio spectra 
with a function that provides the flux density and the 
frequency of the spectral peak, but without any physical 
content. Since most of the sources are optically thin at 
43 GHz, the new flux density measurements at such 
frequency, which was not available in the previous epochs, 
provide very tight constraints on the fits and a better 
determination of the peak. 

In Fig. [TJ we show the radio spectra of all the sources 
observed at the various epochs. 

Following the approach by Torniaincn et al. (2005) we 
compute the spectral indices abciow and a a bovc of the 
overall spectrum (Col. 14 and 15 of Table \2§ fitting a 
straight line to the parts below and above the spectral 
peak respectively. We consider "flat" those sources with 
both abciow > —0.5 and a < 0.5 (where S cx v~ a ). In a 
few sources, depending on the peak frequency, we could 
fit either aboiow or a a bove , in order to avoid the flattening 
near the spectral peak. In this case, sources with the 
spectral index in the range of -0.5 and 0.5 are considered 
flat-spectrum objects. We find that 18 objects, labeled 
with an "F" in Column 9 of Table El no longer show 
the convex spectrum. Such sources are labelled "flat" in 
Column 5 and 6 of Table [31 and are definitely classified 
as "blazar" , and removed from the sample of genuine 



young radio sources. Six of these sources (J0217+0144, 
J0329+3510, J0357+2319, J2123+0535, J2320+0513 and 
J2330+3348) were already found with a flat spectrum in 
the second epoch observations, and already rejected from 
the sample of candidate HFPs (Tinti et al. !2005p . 



3.2. Variability index 

As already mentioned in Section 3.1, given the large 
Doppler factors characterizing the blazar jets, the flux 
density variability can substantially modify the observed 
spectrum of blazars on short timescales. 
On the contrary, young radio sources should not display 
any significant flux density variability and they can be 
considered as the least variable class of extragalactic 
radio sources (O'Dea 1998]), with a mean variation of ~ 
5% (Stanghellini et al. I2005|) . Therefore, genuine young 
HFP objects, considered to be newly born radio sources, 
should not display significant variability. 
We analyze the variability of the sources in terms of the 
quantity: 

TO ^ of 

which is a multi-epoch generalization of the variability in- 
dex defined by Tinti et al. (|2005p . Si is the flux density 
at the i-th frequency measured at one epoch, while S% is 
the mean value computed averaging the flux density at 
the i-th frequency measured at all the available epochs; <ii 
is the error on Si — Si, and to is the number of sampled 
frequencies. 

Columns 7 and 8 of Table [3] report the variability between 
each new epoch and the mean value obtained by averaging 
all the available epochs. The variability index V has been 
computed for each single new epoch, rather than consid- 
ering all the epochs together. In this case the availability 
of two distinct values better indicates the presence of flux- 
density bursts. 

Comparing the variability distributions computed from 

Eq.Q]for all the sources of each epoch, the KS test does not 

detect any significant difference (>99%). This means that 

the flux density of the majority of the observed sources 

has not changed its behaviour with time. 

From the comparison of the multi-epoch radio spectra we 

find: 

— 18 sources with a convex spectrum at the first epoch 
(Dallacasa et al. I2000|) show a flat spectrum in at least 
one of the subsequent epochs. They have an "F" in 
Column 9 of Table [3] and they are rejected from the 
sample; 

— Of the aforementioned sources the quasar J2320+0513 
shows a continuous alternation between a flaring-phase 
with a convex-shape spectrum, and a quiescent-phase 
where the spectrum is flat. In particular, during the 
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Source 


Id. 


z 


code 


Si. 4 


Si. 7 


S4.5 


S4.9 


Ss.i 


S8.4 


S15.3 


S22.2 


S43.2 


Q^bclow 


C^abovc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


J0003+2129 


G 


0.455 


d 


100 


125 


258 


262 


233 


227 


146 


70 


- 


-0.8 


0.9 








e 


84 


- 


251 


255 


226 


220 


134 


70 


16 


-0.9 


1.3 


J0005+0524 


Q 


1.887 


c 


166 


174 


212 


206 


169 


165 


122 


84 


- 


- 


0.6 








e 


168 


- 


213 


207 


170 


168 


129 


90 


46 


- 


0.6 


J0037+0808 


G 




c 


98 


119 


287 


290 


273 


267 


197 


139 


- 


-0.9 


0.5 








e 


95 


- 


287 


293 


278 


274 


210 


148 


79 


-0.9 


0.6 


J0111+3906 


G 


0.668 


c 


509 


597 


1319 


1301 


965 


918 


449 


273 


111 


-0.8 


1.1 


J0116+2422 


EF 




c 


106 


122 


245 


244 


240 


226 


177 


121 


- 


-0.7 


0.4 








e 


99 


- 


258 


238 


240 


225 


177 


129 


- 


- 


0.4 


J0217+0144 


Q 


1.715 


c 


577 


514 


337 


326 


268 


263 


244 


239 


- 


- 


0.3 








e 


571 


- 


326 


315 


257 


264 


256 


254 


- 


- 


0.3 


J0329+3510 


Q 


0.5 


d 


410 


443 


603 


611 


586 


578 


523 


447 


- 


-0.3 


0.2 








g 


563 


- 


644 


625 


540 


535 


468 


401 


302 


- 


0.2 


J0357+2319 


Q 




d 


116 


124 


242 


255 


309 


312 


382 


373 


- 


-0.5 


- 








g 


- 


- 


185 


187 


208 


210 


239 


232 


212 


-0.2 


0.1 


J0428+3259 


G 


0.479 


d 


148 


195 


513 


524 


539 


531 


421 


269 


- 


-0.8 


0.6 








g 


172 


199 


497 


512 


522 


520 


394 


253 


104 


-0.7 


0.9 


J0519+0848 


EF 




c 


269 


263 


394 


393 


373 


372 


365 


336 


- 


-0.4 


0.1 








g 


262 


273 


314 


320 


322 


321 


372 


382 


348 


-0.1 


0.1 


J0625+4440 


BL 




c 


173 


186 


223 


222 


219 


218 


197 


174 


- 


-0.2 


0.1 








g 


121 


131 


180 


182 


187 


187 


180 


161 


121 


-0.2 


0.2 


J0638+5933 


EF 




c 


249 


315 


605 


619 


680 


678 


701 


653 


- 


-0.5 


- 








g 


300 


325 


597 


608 


649 


651 


706 


643 


412 


-0.4 


0.5 


J0642+6758 


Q 


3.180 


c 


226 


272 


436 


429 


343 


332 


205 


131 


- 


-0.5 


0.7 








g 


190 


238 


405 


402 


320 


312 


204 


123 


44 


-0.6 


0.9 


J0646+4451 


Q 


3.396 


c 


444 


660 


2860 


3071 


4068 


4124 


3944 


3184 


- 


-1.3 


0.3 








g 


586 


717 


2894 


3103 


4064 


4094 


4039 


3175 


1996 


-1.1 


0.5 


J0650+6001 


Q 


0.455 


c 


507 


633 


1161 


1150 


994 


975 


674 


467 


- 


-0.7 


0.5 








g 


480 


626 


1150 


1106 


958 


935 


798 


466 


204 


-0.7 


0.7 


J0655+4100 


G 


0.02156 


c 


209 


232 


323 


330 


354 


354 


319 


268 


- 


-0.3 


0.3 








g 


198 


233 


351 


341 


369 


369 


360 


307 


226 


-0.3 


0.3 


J0722+3722 


Q 


1.63 


d 


148 


180 


199 


199 


178 


174 


138 


99 


- 


- 


0.4 








g 


171 


188 


203 


203 


176 


172 


125 


83 


44 


- 


0.7 


J0927+3902 


Q 


0.6948 


d 


2810 


3615 


9545 


9760 


10027 


9937 


8813 


7237 


- 


-0.8 


0.3 


J1016+0513 


Q 




d 


633 


711 


489 


478 


428 


420 


388 


327 


- 


- 


0.3 


J1045+0624 


Q 


1.507 


d 


185 


245 


296 


289 


268 


266 


226 


152 


- 


- 


0.4 


J1148+5254 


Q 


1.632 


a 


108 


- 


396 


414 


460 


450 


411 


277 


- 


-0.9 


0.4 


J1335+4542 


Q 


2.449 


a 


267 


- 


821 


821 


666 


646 


392 


263 


- 


-0.9 


0.7 


J1335+5844 


EF 




a 


299 


- 


745 


744 


727 


726 


585 


449 


- 


-0.8 


0.3 


J1407+2827 


G 


0.0769 


b 


865 


1133 


2519 


2532 


2147 


2071 


1027 


542 


- 


-0.8 


1.0 


J1412+1334 


EF 




b 


191 


248 


346 


341 


286 


277 


189 


121 


- 


-0.5 


0.6 


J1424+2256 


Q 


3.626 


b 


371 


480 


652 


637 


- 


394 


247 


144 


- 


-0.4 


0.9 


J1430+1043 


Q 


1.710 


b 


321 


423 


865 


861 


780 


767 


582 


473 


- 


-0.8 


0.4 


J1505+0326 


Q 


0.411 


b 


382 


428 


608 


620 


620 


610 


515 


468 


- 


-0.3 


0.2 


J1511+0518 


G 


0.084 


b 


92 


123 


569 


608 


801 


811 


763 


573 


- 


-1.3 


0.3 


J1526+6650 


Q 


3.02 


d 


109 


- 


426 


440 


380 


369 


211 


104 


- 


-1.1 


1.0 


J1623+6624 


G 


0.203 


d 


154 




282 


281 


257 


254 


191 


123 




-0.5 


0.5 








i 


162 


179 


265 


263 


227 


221 


164 


122 


80 


-0.4 


0.5 


J1645+6330 


Q 


2.379 


d 


280 




419 


430 


464 


463 


410 


273 




-0.3 


0.5 








i 


270 


288 


391 


399 


430 


432 


401 


322 


192 


-0.3 


0.5 


J1717+1917 


Q 


1.81 


h 


211 


211 


232 


232 


210 


208 


174 


144 


124 


-0.1 


0.3 


J 1735+5049 


G 




h 


445 




925 


943 


898 


888 


678 


485 


262 


-0.6 


0.6 








i 


436 


479 


935 


945 


886 


873 


623 


408 


140 


-0.6 


0.9 


J1751+0939 


BL 


0.322 


h 


980 


1120 


2714 


2848 


3825 


3897 


4556 


4583 


4642 


-0.5 




J1800+3848 


Q 


2.092 


f 


262 


326 


793 


836 


1109 


1128 


1249 


1147 


758 


-0.7 


0.5 



Table 2. The VLA flux density of the 51 candidates HFPs observed with the VLA during the observing runs presented 
in this naner. Col. 1: source name (.72000): Col. 2: ontioal identification: Col. 3: redshift: Col. 4: the observation code. 
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Source 


Id. 


z 


code 


Si. 4 


Si.r 


S4.5 


S4.9 


S8.1 


S8.4 


S15.3 


S22.2 


S43.2 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


J 1840+3900 


Q 


3.095 


f 


123 


147 


161 


163 


165 


165 


165 


164 


140 


-0.1 


0.1 


J1850+2825 


Q 


2.560 


f 


235 


280 


1097 


1185 


1515 


1520 


1402 


1132 


605 


-1.1 


0.6 


J1855+3742 


G 




f 


180 


181 


360 


344 


215 


206 


123 


86 


50 


-0.6 


0.9 


J2021+0515 


Q 




f 


360 


442 


520 


506 


407 


397 


285 


206 


100 


-0.3 


0.7 


J2024+1718 


Q 


1.05 


f 


307 


324 


594 


609 


623 


619 


519 


403 


242 


-0.5 


0.6 


J2101+0341 


Q 


1.013 


c 


483 


493 


497 


498 


552 


555 


704 


724 


- 


-0.2 


- 








f 


431 


- 


498 


508 


665 


687 


954 


1031 


883 


-0.3 


0.2 


J2114+2832 


Q 


2.345 


d 


414 


493 


612 


599 


545 


535 


458 


357 


- 


-0.3 


0.3 


J2123+0535 


Q 


1.878 


f 


2185 


2260 


2831 


2879 


3054 


3057 


3021 


2740 


1992 


-0.1 


0.2 


J2136+0041 


Q 


1.932 


c 


4234 


5231 


10284 


10225 


8763 


8601 


6252 


4401 


- 


-0.7 


0.5 








f 


3752 


4823 


10193 


10150 


8779 


8599 


6343 


4720 


2578 


-0.8 


0.6 


J2203+1007 


G 




c 


107 


156 


315 


311 


240 


231 


125 


67 


21 


-0.9 


1.2 


J2207+1652 


Q 


1.64 


d 


177 


206 


220 


224 


227 


226 


214 


177 




-0.1 


0.2 








c 


268 




242 


248 


246 


246 


218 


183 


151 




0.2 


J2212+2355 


Q 


1.125 


d 


490 


576 


659 


651 


629 


625 


631 


547 




-0.2 


0.1 








c 


507 




634 


644 


682 


684 


695 


624 


517 


-0.1 


0.3 


J2257+0243 


Q 


2.081 


c 


199 


199 


280 


295 


400 


408 


480 


426 




-0.4 


0.3 








c 


180 




291 


306 


420 


432 


517 


466 


330 


-0.5 


0.5 


J2320+0513 


Q 


0.622 


c 


526 


541 


968 


1006 


1144 


1148 


1114 


936 




-0.5 


0.2 








c 


615 




1060 


1076 


1098 


1094 


996 


844 


616 


-0.4 


0.3 


J2330+3348 


Q 


1.809 


d 


274 




498 


511 


596 


603 


620 


525 




-0.3 


0.4 



Table 2. Continued 



second-epoch observations, the source was character- 
ized by a flat spectrum, while during our last observ- 
ing run its spectrum becomes convex again and with 
almost the same flux density as the first epoch. This ex- 
ample points out how important a multi-epoch, multi- 
frequency flux density monitoring is in order to reveal 
blazar objects. 

— 14 sources maintain a convex spectrum at the various 
observing epochs, although with significant flux den- 
sity variability (V >3), and they have a "V" in Column 
9 of Table ED 

— 20 sources preserve the convex spectrum and do not 
show significant flux density variability (V <3), and 
they have an "H" in Column 9 of Table El 



When we compare the variability properties between 
sources with different optical identification by means 
of a KS test, we find that there is a difference (>90%) 
between the variability of galaxies and quasars, as also 
found by Torniainen et al. (|2007[) . supporting the idea 
that radio sources with different optical identification 
represent different radio source populations (Stanghcllini 
et al. I2005| Orienti et al. I2006|) . This difference becomes 
stronger (>99%) if in the KS test we consider sources 
with or without a CSO-like pc-scale morphology (Orienti 
et al. 1200511 . 



3.3. Peak frequency 

So far the anti-correlation (O'Dea & Baum |T997| between 
the peak frequency and the projected linear size has been 
explained mainly in terms of Synchrotron Self- Absorption 
(SSA): as the radio source expands the turnover moves 
progressively to lower frequencies as the result of a 
decreased energy density within the emitting region. 
In Table [3] we report the peak frequency measured at the 
various epochs. A KS test considering all the observed 
sources does not detect any significant (>99%) difference 
among the distributions of the peak frequency at the 
various epochs. This result is expected since the time 
elapsed between the observing runs is too short to detect 
any modification in the spectra of the growing sources. 
However, if we consider individual objects we find that 
most of them have a smaller v v at the subsequent epochs, 
consistent with the evolution models. The median value of 
the peak frequency of the whole sample has continuously 
decreased: v v = 6.7 GHz at the first epoch (Dallacasa et 
al. [2000), v v — 6.3 GHz at the second epoch (Tinti et al. 
2005), and v v = 6.0 GHz at the epochs presented here. 
If we consider young HFP candidates and blazar objects 
separately, we find that the former show a decreasing 
peak frequency, from a median value of 6.0 GHz during 
the first epoch (Dallacasa et al. 2000) to 5.5 GHz at 
the subsequent epochs. The latter have a median peak 
frequency which does not follow a monotonic trend: v p = 
7.4 GHz at the first epoch, v v — 8.6 GHz at the second 
epoch, v v = 6.2 at the third epoch and v v — 6.8 GHz at 
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Fig. 2. Canonical radio spectra computed for radio 
sources with V<3 (top panel); radio sources with V>3 
(bottom panel). The continuous line indicates the least- 
square fit to the data. 



the last epoch. A KS test does find a difference (>99%) 
between the peak frequency distributions between young 
HFP candidates and blazars, although at the first two 
epochs only. 

On the other hand, among the 14 sources with V>3 
(Section 3.2) we find 3 objects in which the turnover 
frequency shows a remarkable change, although the 
overall radio spectrum maintains a convex shape. Two of 
these sources (the quasars J1645+6330 and J2024+1718) 
show a large drop in the peak frequency. We estimate the 
size of the emitting region by two independent methods. 
In one case, we use the relationship from O'Dea (1998): 

logi/p = -0.21 - omiogLS 

which relates the intrinsic spectral peak and the linear 
size. In the other one, we assume synchrotron self- 
absorption theory (Kellerman &; Pauliny-Toth I1981|) , in 



which: 



= B*S$v p 4 (l + z)- 



(2) 



where B is the magnetic field, S p the peak flux density, 
i/p the turnover frequency and z the redshift. For the 
magnetic field we consider the values reported by Orienti 
et al. (|2006j) . obtained assuming equipartition conditions. 
As will be discussed in more detail in the following 
section, if we assume that the magnetic field is frozen 
within an adiabatically-expanding homogeneous region, 
its value can be considered constant during the 5 years 
that elapsed between the most distantly separate observ- 
ing runs. 

The increment in size obtained in both ways are in good 
agreement, and corresponds to an expansion velocity > c, 
which is clearly unrealistic for unboosted young objects. 
We conclude that such sources are beamed objects and 
we reject them as genuine candidate HFPs. 
The other object, the BL Lac J1751+0939, shows an 
increment of the peak frequency (from 8.5 to 29 GHz), 
which can be interpreted in terms of different knots in the 
jet base dominating the radio emission at the two epochs. 
This source has a well known history of flux density and 
spectral variability, and therefore there is no doubt about 
its blazar nature. 

In general, the observed peak frequency of galaxies and 
quasars are similar as a consequence of the selection 
criteria. However, if we consider the intrinsic turnover 
frequency, the KS test finds a clear difference (>99%) 
concerning: 

a) galaxies and quasars; 

b) the parsec-scale structure depending on its having 
the CSO-like morphology typical of young radio sources 
(Orienti et al. l2005|) . 

While in the former case the different distribution is 
easily explained in terms of redshift (quasars are found 
at higher redshifts than galaxies), in the latter, such a 
segregation is likely to be indicating two different source 
populations. 

In order to analyze the spectral shape of our sample 
sources, we construct a canonical radio spectrum for 
sources with either V<3 (Fig. Hi) or V>3 (Fig. [2b). 
Following the work of de Vries et al. (|1997p . the canonical 
radio spectra have been obtained by normalizing the 
observed frequencies and fluxes by the source peak 
frequency and the peak flux density, computed averaging 
all the epochs: 



_ \ " 



1 Jz> c 



(3) 



(4) 



pi 



where v and S are the average normalized frequency 
and flux density respectively, v \ and v p i are the observed 
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frequency and the observed turnover frequency at the 
i-th epoch, while Si and S p j are the observed flux density 
and the observed peak flux density at the «-th epoch, and 
m the epochs available. 

Sources with V<3 or V>3 display different canonical 

spectral shape. The former have a quite convex spectrum, 

with a narrow width (FWHM ~1.7). The spectral indices 

are: athick ~-0.9±0.1 and ccthin ~ 0.7±0.1. 

The latter have a flatter spectral shape, with a FWHM 

~3.6 and spectral indices: athick ~-0.4±0.1 and athin ~ 

0.4±0.1. 



4. Discussion 

A multi-epoch monitoring of the spectral characteristics 
provides information on the evolution of the radio spectra. 
From the analysis of the spectral shape we find 12 sources 
(see Section 3.1) showing a fiat spectrum during at least 
one of the observing epochs presented in this paper, in 
addition to the 7 flat-spectrum objects already rejected 
by Tinti et al. (|2005|) . By comparing the frequency break 
obtained at different epochs we also find that in 3 sources 
(see Section 3.3) the break has noticeably shifted. Such a 
change can be explained only assuming a beamed nature 
for these sources, and therefore we reject them from the 
sample of genuine young HFPs. If we also consider the 7 
sources (5 of them even with a flat spectrum) found with 
a Core- Jet morphology by Orienti et al. (12006|) and then 
rejected, we obtain that only 31 of the 55 sources of the 
sample can still be considered HFP candidates, and they 
are marked in boldface in Table [3l 

Although the identification of contaminant objects is 
the main goal, a multi-epoch monitoring program can 
provide important information on the spectral evolution 
in young radio sources. 

We assume that young radio sources are described by 
a continuous injection model, where the radio emission 
is continuously replenished by a constant flow of fresh 
relativistic particles with a power-law energy distribution. 
However, the overall spectra of such sources, charac- 
terized by a convex shape, show significant departure 
from the expected classical power-law (S oc v~ a ). The 
deviations are well explained by two different phenomena: 
synchrotron self-absorption causes the rising spectrum (S 
oc v^l 2 ) at frequencies below the peak, while synchrotron 
losses steepen the spectrum (S oc i/~("+ - 5 )) at high fre- 
quencies (Pacholczyk 1970). In very young radio sources 
such as the HFPs, such a steepening of the spectrum 
occurs at very high frequency implying that synchrotron 
losses cannot be investigated with our VLA observations. 
In the following discussion we investigate the nature of 
the radio spectrum and its evolution. 



4.1. Adiabatic expansion 

At frequencies below the turnover, the adiabatic expan- 
sion is the main mechanism that influences the spectral 
evolution. 

In this regime we have: 



S(v) oc 



(5) 



where B is the magnetic field, v the frequency and 9 
the angular size of the emitting region. 
As the radio source adiabatically expands, the opacity 
decreases, the turnover moves to lower frequency and the 
flux density increases. 

Is it possible to detect such an increment in our multi- 
epoch spectra? 

First, we assume that the radio emission is due to a 
homogeneous component that is adiabatically expanding 
at a constant rate: 



t + At 
to 



(6) 



where 9o is the angular size at the time to, and 9 at the 
time to+At, and the magnetic field is frozen in the plasma: 



B = B n 



to 



t a + At 



(7) 



where Bo is the magnetic field at the epoch to, while 
B is at the time to + At. 

From Equation the flux density at the frequency v at 
a given time t is: 



The turnover frequency is: 



v t ,o oc B E 



(8) 



(9) 



where Eq is the energy of the relativistic particles. During 
the expansion of the homogeneous synchrotron-emitting 
region E oc t -1 (Pacholczyk 1 1970")) . which implies that 



"t,i = ft,o 



to 



t + 



At 



(10) 



If t — to=At <C! to, the flux density and the turnover 
frequency can be considered approximately constant. We 
consider At ~ 5 years (the time elapsed between the first 
and last observing run), and to > 15 years, since all the 
sources are part of the 87GB sample. If to = 20 years, 
we expect a flux density increment of a factor of ~ 2 in 
the optically thick part of the spectra. Such an increment 
becomes undetectable for sources with ages of ~ 100 
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Source 


Morph. 


Vepl 


Vep2 


^ep3 


fep4 


Vepl 


V cp2 


Var. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


J0003+2129 


cso 


5.7±0.1 


5.4±0.1 


5.2±0.9 


5.29±0.02 


2.67d 


3.32e 


H 


J0005+0524 


cso 


4.13±0.09 


3.40±0.09 


3.6±0.7 


3.5±0.7 


1.65c 


2.1e6 


H 


J0037+1109 


cso 


5.9±0.1 


6.2±0.2 


6.2±0.7 


6.0±0.7 


1.03c 


1.47e 


H 


J0111+3906 


cso 


4.76±0.06 


4.68±0.07 


4.66±0.19 




4.63c 




H 


J0116+2422 


Un 


5.1±0.1 


6.3±0.3 


6.0±0.8 


6.0±1.1 


4.93c 


3.73e 


V 


J0217+0144 


Un 


18±1 


flat 


flat 


flat 


94.72c 


94.10e 


F 


J0329+3510 


CJ 


6.7±0.3 


flat 


flat 


flat 


4.31d 


24.43g 


F 


J0357+2319 


Un 


12±1 


flat 


flat 


flat 


12.67d 


42.48g 


F 


J0428+3259 


CSO 


7.3±0.2 


6.8±0.2 


6.9±0.5 


6.55±0.38 


1.12d 


2.40g 


H 


J0519+0848 


Un 


>22 


7.4±0.6 


7.26±0.01 


flat 


13.28c 


19.84g 


F 


J0625+4440 


Un 


13±2 


7.4±1.0 


flat 


flat 


18.18c 


62.38g 


F 


J0638+5933 


CSO 


12±2 


9.2±0.7 


11.35±0.30 


10.22±0.20 


1.79c 


3.74g 


H 


J0642+6758 


Un 


4.5±0.1 


4.08±0.08 


4.41±0.46 


4.41±0.54 


0.71c 


9.18g 


V 


J0646+4451 


MR 


15±2 


11.2±0.6 


10.36±0.09 


10.65±0.05 


19.19c 


32.44g 


V 


J0650+6001 


CSO 


7.6±0.3 


5.2±0.2 


5.23±0.18 


5.45±0.15 


2.89c 


3.12g 


V 


J0655+4100 


Un 


7.8 




flat 


flat 


0.46c 


3.73g 


F 


J0722+3722 


MR 


4.3 




4.0±0.7 


3.4±0.6 


2.67d 


4.76g 


H 


J0927+3902 


CJ 


6.9 




8.3±0.1 




15.45d 




V 


J1016+0513 


Un 


7.1 




flat 




72.80d 




F 


J1045+0624 


Un 


3.7 




4.7±0.5 




5.75d 




H 


J1148+5254 


CSO 


8.7 




7.9±0.7 




9.74a 




V 


J1335+4542 


cso 


5.1±0.1 


4.9±0.1 


5.1±0.4 




4.64a 




H 


J1335+5844 


cso 


6.0±0.2 


5.5±0.1 


6.50±0.34 




2.84a 




H 


J1407+2827 


cso 


5.34±0.05 


5.01±0.08 


4.95±0.01 




2.02b 




H 


J1412+1334 


Un 


4.7±0.1 


4.18±0.09 


4.3±0.5 




0.97b 




H 


J1424+2256 


Un 


4.13±0.07 


3.94±0.06 


3.7±0.3 




8.80b 




V 


J1430+1043 


MR 


6.5±0.2 


5.7±0.1 


6.2±0.1 




1.37b 




H 


J1505+0326 


Un 


7.1±0.4 


6.8±0.4 


flat 




30.13b 




F 


J1511+0518 


CSO 


11.1±0.4 


10.8±0.4 


10.10±0.04 




10.24b 




V 


J1526+6650 


MR 


5.7±0.1 


5.5±0.1 


5.5±0.9 




2.04d 




H 


J1623+6624 


Un 


6.0±0.2 


6.0±0.2 


5.1±0.8 


4.5±0.5 


0.87d 


11.161 


V 


J1645+6330 


Un 


14±2 


10.1±0.7 


6.0±0.1 


6.8±0.3 


12.38d 


21.611 


V 


J1717+1917 


Un 


11.5 




flat 




81.23h 




F 


J1735+5049 


CSO 


6.4±0.2 


6.3±0.3 


5.6±0.2 


5.6±0.2 


0.27h 


2.07i 


H 


J 175 1+0939 


CJ 


8.5 




29.02±0.04 




8.75h 




V 


J1800+3848 


Un 


17±3 


13±1 


13.1±0.1 




1.92f 




H 


J1840+3900 


Un 


5.7±0.5 


5.2±0.4 


flat 




6.87f 




F 


J1850+2825 


MR 


9.1±0.3 


9.5±0.3 


9.9±0.4 




2.15f 




V 


J1855+3742 


CSO 


4.00±0.07 


3.81±0.06 


4.0±0.5 




1.32f 




H 


J2021+0515 


CJ 


3.75±0.08 


4.5±0.1 


3.7±0.3 




9.94f 




V 


J2024+1718 


Un 


14±2 


8.6±0.4 


6.6±0.4 




14.51f 




V 


J2101+0341 


Un 


17±2 


3.7±0.2 


flat 


flat 


28.29c 


17.00f 


F 


J2114+2832 


CJ 


9.8 




flat 




34.76d 




V 


J2123+0535 


CJ 


18±4 


flat 


flat 




112. 74f 




F 



Table 3. Peak frequencies and flux density variability between the different epochs of the 51 sources observed with the 
VLA during the observing runs presented in this paper. The variability indices have been computed by comparing the 
flux densities measured at each epoch with the mean value obtained by averaging all the available observing epochs, 
as described in Section 3.2 (see Col. 7 and 8). Col. 1: source name (J2000); Col. 2: Morphological classification from 
pc-scale information (Orienti et al. 120061) Col. 3 and 4: peak frequency of the first two epochs (Dallacasa et al. 120001 
Tinti et al. 2005); Col. 5 and 6: peak frequency from these new VLA data. Col. 7 and 8: the variability index computed 
by comparing the flux densities measured in 2003/2004 and presented in this paper with those averaged on all the 
available epochs, as described in Section 3.2. Since the data flux densities used in Column 7 and 8 were collected 
during different observing runs, near the value of the variability we report also the observation code from Table [TJ 
which indicates for each source the date of the observations. Col. 9: The classification of the source spectra: H= sources 
without significant variability, see Section 3.2; V= variable (i.e. with a variability index V>3, as defined in Section 3.2 
at least during one epoch); F= flat spectrum (i.e. sources with a flat spectrum at least during one epoch; see Section 
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Source 


Morph. 




^ep2 


^ep3 


^cp4 


Vepl 


V C p2 


Var. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


J2136+0041 


CJ 


5.0 




5.52±0.01 


5.73±0.02 


2.88c 


0.37f 


H 


J2203+1007 


cso 


4.86±0.07 


5.0±0.1 


4.6±0.7 




1.94c 




H 


J2207+1652 


Un 


7.4±0.3 


3.5±0.3 


flat 


flat 


36.01d 


35.45e 


F 


J2212+2355 


Un 


13±2 


9±1 


flat 


flat 


43.10d 


40.29e 


F 


J2257+0243 


Un 


>22 


>22 


flat 


13.6±0.3 


3.27c 


1.92e 


F 


J2320+0513 


Un 


5.4±0.2 


fiat 


9.2±0.2 


flat 


23.84c 


5.75e 


F 


J2330+3348 


MR 


5.6±0.3 


flat 


flat 




18.84d 




F 



Table 3. Continued. 



years. On the other hand, the v t should have decreased 
by about 70% of its previous value. If to ~ 100 years, the 
decrement should be less than 20% of its previous value, 
i.e. a turnover frequency of 22 GHz moves to 18 GHz, 
making the detection rather difficult due to the poor 
frequency coverage available. This result is based on the 
strong assumptions that the magnetic field is frozen in a 
homogeneous component which is adiabatically expand- 
ing. However, it is possible that either the time-variable 
magnetic field is not frozen in the plasma, or the radio 
emitting region is not homogeneous. In both cases the 
flux density increment and the frequency decrement are 
expected to be even less significant. 

Therefore, we conclude that adiabatic expansion is not 
able to produce any significant changes in the optically 
thick part of the spectrum detectable in our multi-epoch 
observations, perhaps with the exception of a few sources, 
such as J0650+6001 (Section 4.3) where a tentative flux 
density increment at a frequency below the spectral peak 
has been found, or the "faint" HFP J1459+3337 (Orienti 
& Dallacasa, submitted). 



have different turnover frequencies, and, thus, different 
parts of the overall spectra are influenced more by one 
component, rather than the other, as is easily seen in the 
case of J1511+0518. As a consequence, changes in the 
total spectra may strongly depend on the evolution of the 
emission of each single component. 

The variation of the total spectrum that occurred over 5 

years makes J0650+6001 an interesting case study. 

If we compare the overall spectra measured at different 

epochs, we find that at 1.4 GHz, well below the peak 

frequency, the flux density appears to be increasing, 

although within the errors, as expected by adiabatic 

expansion. 

On the other hand, at frequencies higher than the peak 
the flux density decreases. During the first epoch of 
observation, the optically thin spectrum was described by 
a power-law with a ~ 0.3, while during the subsequent 
observing runs it displayed a ~ 0.7. Such a steepening 
cannot be directly related to a shift of the break frequency 
of the overall spectrum, but is rather due to the evolution 
of the spectra of the individual components. 



4.2. The nature of the overall spectra 

In the previous Sections we have assumed that the radio 
emission is due to a partially opaque electron synchrotron 
radiation which originates within a homogeneous region. 
Variations in the opacity throughout the emitting region 
lead to an overall spectrum which is the superposition of 
the spectra of many single regions. 

Figure [3] shows the overall spectrum (diamonds) of 
J0650+6001 and J1511+0518, together with the spectra 
of each single component (squares and crosses), as 
revealed by VLBI observations (Orienti et al. I2006|) . In 
both sources, pc-scale resolution images show that the 
radio emission originates within two main components, 
which can be interpreted as the lobes of a mini radio 
source. 

The overall spectra are easily explained as the result of 
the superposition of the spectra of individual compo- 
nents. In both sources, the two radio emitting regions 



5. Conclusions 

We have presented the results of new epochs of high- 
sensitivity simultaneous multi-frequency VLA obser- 
vations of a sample of young HFP candidates. By 
considering the spectral shape and the change of the 
frequency break we find 15 sources which turned out to 
be blazar objects, in addition to the 7 already found with 
a flat spectrum by Tinti et al. ( 2005|) . If we consider also 
the 7 sources rejected by Orienti et al. (I2006|) on the 
basis of their Core- Jet morphology, and 5 of them also 
with a flat spectrum, we find that 24 of the 55 sources 
from the HFP sample are contaminant blazar objects, 
and thus only 31 objects (~56%) can still be considered 
HFP candidates. 

The comparison of the variability properties with the 
optical identification have shown that quasars and galax- 
ies display different characteristics. If in the case of the 
turnover frequency distribution such a segregation is due 
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Fig. 3. The spectra of the sources J0650+6001 (top) 
and J1511+0518 (bottom): diamonds represent the over- 
all spectrum obtained with these new VLA observations, 
while crosses and squares represent the spectra of the two 
source components (Xiang et al. I2002I Xu et al. I1995I 
Orienti et al. f2006|) . In the case of J0650+6001, triangles 
show the overall spectrum from Dallacasa et al. (2000). 
The asterisks represent the sum of the spectra of the sin- 
gle components. In the case of J1511+0518 not all the flux 
density could be recovered by VLBI images (Orienti et al. 
[20061) . 



to redshift ranges (quasars are found at systematically 
higher redshifts than galaxies), the difference in the 
variability index likely reflects two kinds of radio source 
populations. An even stronger segregation in the flux- 
density and spectral-shape variability is found between 
objects with or without a CSO-like morphology. 
Since polarization properties are another effective tool 
to discriminate young radio sources from contaminant 
blazars, polarimetric data concerning the sources of the 
"bright" HFP sample have been studied, and they will 
be presented in a companion paper (Orienti & Dallacasa 
2007, Paper II). 



Making use of the several epochs available, we have tried 
to investigate the spectral evolution of the radio emission 
of genuine young radio sources. 

At frequencies below the peak the spectral evolution 
is dominated by adiabatic expansion, and we would 
expect a flux density increment, while at frequencies 
above the spectral peak, in addition to the decrease of 
the flux density due to adiabatic expansion, synchrotron 
losses also may play a role. Considering the time elapsed 
between the two farthest spaced observing runs (5 years), 
we conclude that a variation in the optically-thick part 
of the overall spectrum would have been detected only if 
the radio emission would have started very recently, i.e. 
~ 20 years, which is statistically unlikely for the sources 
in the sample. 

On the other hand, as seen in the analysis of the spectra 
of individual objects, there are a few cases, such as the 
source J0650+6001, with spectral variations. However, 
if we consider that the radio emission comes from a 
single homogeneous region, we are not able to explain 
such changes. Combining the low-resolution, flux density 
information provided by VLA observations with the 
pc-scale resolution from VLBA data, we can describe the 
total radio spectrum as the result of the superposition of 
the spectra of different emitting regions. 
In general, to determine the evolution of the radio 
emission, we must resolve each single source component. 
For this, new simultaneous multi-frequency VLBA ob- 
servations in both the optically thick and thin parts of 
the spectrum have been carried out for 5 genuine HFPs, 
with a CSO-like morphology. The analysis of the radio 
spectrum of each single component will enable us to set 
strong constraints on the fate of the radio emission. 

Acknowledgements. We thank the referee Merja Tornikoski 
for carefully reading the manuscript and valuable sugges- 
tions. The VLA is operated by the U.S. National Radio 
Astronomy Observatory which is a facility of the National 
Science Foundation operated under a cooperative agreement 
by Associated Universities, Inc. This work has made use 
of the NASA/IPAC Extragalactic Database (NED), which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 

References 

Alexander, P. 2000, MNRAS, 319, 8 

Bicknell, G.V., Dopita, M.A., O'Dea, C.P.O. 1997, ApJ, 485, 
112 

Browne, I.W.A., Patnaik, A.R., Wilkinson, P.N., Wrobel, J. A. 

1998, MNRAS, 293, 257 
Dallacasa, D., Stanghellini, C, Centonza, M., Fanti, R. 2000, 

A&A, 363, 887 
Dallacasa, D. 2003, PASA, 20, 79 

de Vries, W.H., Barthel, P.D., O'Dea, OP. 1997, A&A, 321, 
105 



12 



M. Orienti et al.: Constraining the nature of High Frequency Peakers 



Fanti, C, Fanti, R., Dallacasa, D., Schilizzi, R.T. et al. 1995, 
A&A, 302, 317 

Kameno, S., Horiuchi, S., Shen, Z.-Q. et al. 2000, PASJ, 52, 
209 

Kellerman, K.I., Pauliny-Toth, I.I.K. 1981, ARA&A, 19, 373 

Marecki, A., Spencer, R.E., Kunert, M. 2003, PASA, 20. 46 

Murgia, M. 2003, PASA, 20, 19 

O'Dea, CP, Baum, S.A. 1997, AJ, 113, 148 

O'Dea, CP 1998, PASP, 110, 493 

Orienti, M., Dallacasa, D., Tinti, S., Stanghellini. C 2006, 
A&A, 450, 959 

Orienti, M., Dallacasa, D., Stanghellini, C 2007, A&A, 461, 
923 

Orienti , M., Dallacasa, D. 2007, A&A, submitted 
Pacholczyk, A.G. 1970, Radio Astrophysics, (San Francisco: 

Freeman & Co.) 
Patnaik, A.R., Browne, I.W.A., Wilkinson, P.N., Wrobel, J.M. 

1992, MNRAS, 254, 655 
Polatidis, A.G., & Conway, J.E. 2003, PASA, 20, 69 
Readhead, A.C.S., Taylor, G.B., Xu, W., Pearson, T.J. et al. 

1996, ApJ, 460, 612 
Snellen, I.A.G., Schilizzi, R.T., Miley, G.K. et al. 2000, 

MNRAS, 319, 445 
Stanghellini, C, O'Dea, CP, Dallacasa, D. et al. 2005, A&A, 

443, 891 

Tinti, S., Dallacasa, D., de Zotti, C, Celotti, A., Stanghellini, 

C. 2005, A&A, 432,31 
Torniainen, I., Tornikoski, M., Terasranta, H., Aller, H.D. 2005, 

A&A, 435, 839 

Torniainen, I., Tornikoski, M., Lahteenmaki, A., et al. 2007, 
A&A, 469, 451 

Wilkinson, P.N., Polatidis, A.G., Readhead, A.C.S., Xu, W., 

Pearson, T.J. 1994, ApJ, 432, 87 
Wilkinson, P.N., Browne, I.W.A., Patnaik, A.R. et al. 1998, 

MNRAS, 300, 790 
Xiang, L., Stanghellini, C, Dallacasa, D., Haiyan, Z. 2002, 

A&A, 385, 768 

Xu,W., Readhead, A.C.S., Pearson, T.J., Polatidis, A.G., 
Wilkinson, P.N. 1995, ApJS, 99, 279 



M. Orienti et al.: Constraining the nature of High Frequency Peakers 



13 




. 1. Continued. 



14 



M. Orienti et al.: Constraining the nature of High Frequency Peakers 




10.0 100.0 




J2207+1652 


Q 




J2212+2355 


Q 


























0.01 







10.0 100.0 



10.0 100.0 1000.0 0.1 



10.0 100.0 1000.0 



***** 



10.0 100.0 1000.0 0.1 



10.0 100.0 1000.0 



Fig. 1. Continued. 



